A quasi-physical family of gravity-wave templates for precessing binaries of spinning 
compact objects: II. Application to double-spin precessing binaries 
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The gravitational waveforms emitted during the adiabatic inspiral of precessing binaries with two 
spinning compact bodies of comparable masses, evaluated within the post-Newtonian approximation, 
can be reproduced rather accurately by the waveforms obtained by setting one of the two spins to 
zero, at least for the purpose of detection by ground-based gravitational- wave interferometers. Here 
we propose to use this quasi- physical family of single-spin templates to search for the signals emitted 
by double-spin precessing binaries, and we find that its signal-matching performance is satisfactory 
for source masses (mi, 1x12) € [3, 15]Mq x [3, 15]Mq. For this mass range, using the LIGO-I design 
sensitivity, we estimate that the number of templates required to yield a minimum match of 0.97 is 
~ 320, 000. We discuss also the accuracy to which the single-spin template family can be used to 
estimate the parameters of the original double-spin precessing binaries. 
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I. INTRODUCTION 

An international network of long-baseline laser- 
interferometric gravitational-wave detectors, consisting 
of the Laser-Interferometer Gravitational-wave Observa- 
tory (LIGOLQ, of VIRGO @, of GEO 600 and of 
TAMA 300 4], has by now begun science operations. 
VIRGO is in its commissioning phase, while LIGO has 
already completed three science runs (SI on August- 
September 2002 H], S2 on February-April 2003, and S3 
on October 2003-January 2004; SI and S3 were in coinci- 
dence with GEO 600) with increasing sensitivity and sta- 
bility. The analysis of SI data has been completed, yield- 
ing new upper limits on event rates for various classes of 
astrophysical sources 0; the data from S2 and S3 are 
still being analyzed. LIGO is expected to reach its full 
design sensitivity in 2005. 

Compact binaries consisting of black holes (BHs) and 
neutron stars (NSs) are among the most promising 7] 
and best understood sources for such gravitational-wave 
(GW) interferometers, which can observe the waves emit- 
ted by the binaries in the adiabatic-inspiral re gim e, where 
post-Newtonian (PN) calculations H H [TiJ [H 111 13 
fiH HH ll6L IttT ] are appropriate to describe the orbital 
dynamics and predict the gravitational waveforms. 

Very little is known about the statistical distributions 
of BH spins in compact binaries: the spins could be large, 
and they need not be aligned with the orbital angular 
momentum. When this is the case, spin-orbit and spin- 
spin interactions can cause the rapid precession of the or- 
bital plane of the binary, and thus significant modulations 
of the emitted GWs, as it was shown by Apostolatos, 
Cutler, Sussman and Thome (ACST) and later by 
Apostolatos [l^. These modulational effects should be 
included in the theoretical waveform models (templates) 
used in matched-filtering GW searches. However, using 



template banks parametrized by all the relevant phys- 
ical parameters (the masses, the spins, the angles that 
describe the relative orientations of detector and binary, 
and the direction of propagation of GWs to the detector) 
would make such searches extremely computationally in- 
tensive. 

One possibility to reduce the computational cost is the 
adoption of smaller detection template families (DTF), 
which capture the essential features of the true wave- 
forms, but depend on a smaller number of parameters. 
A DTF for precessing binaries was first proposed by 
Apostolatos [H H2, building on the analysis of pro- 
cessional dynamics of Refs. [3 113 • However, accord- 
ing to Apostolatos' own estimates and to Grandclement, 
Kalogera and Vecchio's later tests ,22f , the computational 
resources required by the Apostolatos DTF are still pro- 
hibitive, and its signal-matching performance is unsatis- 
factory. The latter is improved in a modified version of 
the DTF [H , which adds 5-like spikes in the wave- 
form phase. 

Buonanno, Chen and Vallisneri 25, henceforth BCV2] 
investigated the dynamics of precession further, and pro- 
posed a new convention to write the dominant quadrupo- 
lar contribution to GW emission. In this convention, the 
oscillatory effects of precession are isolated in the evolu- 
tion of the GW polarization tensors, which are combined 
with the detector's antenna patterns to yield its response. 
As a result, the response can be written as the product 
of a carrier signal and a modulational correction, which 
can be handled using an extension of Apostolatos' treat- 
ment of precessional effects. BCV2 cast these waveforms 
into a mathematical structure that allows searching au- 
tomatically and economically over all the parameters re- 
lated to precessional modulations, except for a single pa- 
rameter that describes the timescale of modulation. The 
BCV2 DTF has reasonable computational requirements 
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and good signal-matching capabilities. However, espe- 
cially for binaries with high, comparable masses, it has 
the shortcoming that a large number of unphysical wave- 
forms are automatically included in GW searches (albeit 
at no extra computational cost), increasing the probabil- 
ity of false alarms triggered by noise. [This shortcoming 
is unfortunately but unavoidably characteristic of DTFs, 
which replace a description in terms of physical source 
parameters by one in terms of phenomenological signal 
parameters.] 

This paper is the second in a series (begun with Ref. 
[2(|) written to investigate the possibility of searching 
for precessing binaries using a physical family of signal 
templates (a PTF) computed from the PN equations of 
motion. Although at first sight the number of physical 
parameters necessary to describe a waveform is large, we 
were able to reduce the effective dimensionality of the 
template family using the insight developed in the con- 
struction of DTFs. As mentioned above, BCV2 [25L Sec. 
VI D] established that it is possible to search easily over 
most of the parameters related to the kinematics of pre- 
cession (such as the orientation of the detector and of 
the binary as a whole, the direction of GW propagation, 
and the initial orbital phase). In effect, these extrinsic 
parameters can be incorporated in the detection statis- 
tic, while single "templates" [2^ remain functions only 
of the masses of the binary components, of the magni- 
tudes of their spins, and of the relative angles between 
the spins and the orbital angular momentum at a fiducial 
freque ncy . Under the assumption of circular adiabatic in- 
spiral 28] , seven such intrinsic parameters are needed for 
a generic binary where both spins are important (hence- 
forth, a double-spin binary); four are needed for a binary 
where only one body has significant spin (henceforth, a 
single-spin binary). (See Sec. [H] for the distinction be- 
tween extrinsic and intrinsic parameters.) 

In Ref. [2(|, we demonstrated the feasibility of a PTF 
search for single-spin binaries: we described a two-stage 
algorithm to search over the extrinsic parameters (the 
first stage emphasizes computational efficiency, but re- 
tains some unphysical waveforms; the second stage, ap- 
plied only to first-stage triggers, restricts the possible 
search outcomes to physical configurations). Using this 
algorithm, we tested a four-parameter PTF for binaries 
with (mi,m 2 ) € [7M Q ,12M Q ] x [1M Q ,3M Q ], where the 
assumption of a single significant spin is justified. We 
found that ~ 76,000 templates are required for a search 
in this mass range, for a minimal match of 0.97 (see Sec. 
[n]for a definition of minimal match); under the assump- 
tion of Gaussian, stationary noise, we also found that the 
detection threshold required for a given false-alarm prob- 
ability is lower in the PTF search than in a DTF search 
with the same number of intrinsic parameters. 

In this paper we examine PTF searches for the more 
general class of double-spin binaries. Although in this 
case we have seven intrinsic parameters, they are not all 
essential in determining the waveforms. This is str ictly 
true in two limits. First, as it was realized by ACST [18|, 



if the two binary masses are equal, and if spin-spin in- 
teractions are ignored, the same orbital evolution can 
be replicated by giving the total spin to one of the ob- 
jects. Indeed, for the mass ranges of interest to ground- 
based interferometers, spin-spin effects contribute mildly 
to the binding energy and to the PN GW flux, even 
close to the last stable orbit. Second, if the mass ra- 
tio i] = m^rail (mi +7712) 2 is very low (as it was assumed 
in Ref. |2g), the spin of the less massive object can be 
ignored. In addition, as investigated by BCV2 (and less 
systematically by Kidder ^J), the dynamics of double- 
spin binaries with generic mass ratios show features sim- 
ilar to those described by ACST for single-spin binaries. 

These arguments have led us to conjecture that single- 
spin waveforms may always be sufficient to approxi- 
mate double-spin waveforms, at least for the purpose of 
GW searches with ground-based interferometers. Since 
the single-spin parameters that best reproduce a given 
double-spin signal might not be in the physical range for 
a true single-spin binary (for instance, because the spin 
of one of two objects must do the work of two, it might 
have to exceed the maximal spin allowed for BHs), the 
single-spin family should be called quasi-physical, but we 
shall continue to use "PTF" loosely. In the rest of this pa- 
per, we present evidence that our conjecture is correct for 
the mass range (mi, m 2 ) € [3M©, 15M©] x [3M Q , 15M Q ], 
and we examine the computational requirements and the 
parameter-estimation performance of a single-spin PTF 
search for double-spin binaries in this mass range. 

This paper is organized as follows. In Sec. II, we pro- 
vide a short glossary for the notions and quantities of 
matched-filtering GW searches (some standard, some de- 
veloped in Ref. |2j|) that are needed later. In Sec. Ill A, 
we review the adiabatic PN dynamics of double-spin bi- 
naries; in Sec. Ill B, we describe our family of quasi- 
physical single- spin templates, and we evaluate their 
signal-matching capabilities against double-spin binaries 
with maximal spins (where precessional effects are ex- 
pected to be strongest); in Sec. Ill C, we study the ro- 
bustness of adiabatic PN waveforms for binaries with 
high, comparable masses; in Sec. Ill D, we discuss some 
features of double-spin binary dynamics that help to ex- 
plain the good signal-matching performance of single- 
spin templates. In Sec. IV, we estimate the number of 
templates required in a single-spin PTF search in our 
mass range of interest. In Sec. V, we investigate the ex- 
traction of the physical parameters of the double-spin 
binary using single-spin templates. Last, in Sec. VI we 
summarize our main conclusions. 



II. A GLOSSARY OF MATCHED-FILTERING 
GW DETECTION 

In this paper we adopt the standard formalism of 
matched-filtering GW detection, as summarized in Ref. 
|17| (which includes an extensive bibliography), and as 
extended in Ref. [2(| to a special treatment of extrin- 
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sic and intrinsic parameters. Here we provide a glossary 
of the notions used in this paper, with pointers to their 
definitions in Refs. [HSl]. 

Templates h(X A ): Theoretical models of GW signals, 
parametrized by one or more template parameters X A . 
A continuous template family {h(X A )} defines a smooth 
submanifold in signal space. 

Noise inner product (f,g): Eq. (1) of Ref. 17}. A 
measure of the closeness of two signal, as given by a 
correlation product weighted by the power spectral den- 
sity of noise; throughout this paper we adopt the LIGO-I 
one-sided noise power spectral density S n given by Eq. 
(28) of Ref. 0| ■ A normalized template h(X A ) has 
(h(X A ),h(X A )) = l. 

Match: Inner product of two normalized signals. The 
mismatch is one minus the match. 

Overlap p(s, h(X A )): Inner product of a signal s with 
the normalized template h(X A ). 

Detection statistic: Figure of merit compared with a 
detection threshold to decide whether the signal modeled 
by a template h(X A ) is present in the detector output o. 
For Gaussian, stationary noise, the overlap pip, h(X A )) 
is the optimal statistic that minimizes the probability of 
false dismissal for a given probability of false alarm (set 
by the detection threshold). In this context p is also 
known as the signal-to-noise ratio (S/N) of the detector 
output after filtering by the template h(X A ). The corre- 
sponding detection statistic for an entire template family 
{h(X A )} is the maximized overlap m&x^A p(o, h(X A )). 

Fitting factor FF 0: Eq. (20) of Ref. [l7|. Match 
between a template in a target family (representing ac- 
tual physical signals) and a template in a search family, 
maximized over all the parameters of the search family. 
The FF (a value between and 1) characterizes the ef- 
fectualncss of the search family in reproducing signals 
modeled by the target family: using an imperfect family 
means that only a fraction FF of the available S/N is 
recovered, reducing the number of true events that pass 
the detection threshold. The maximized match induces 
a (many-to-one) map between the space of target param- 
eters and the space of search parameters (see Sec. V on 
the systematic errors in parameter estimation induced by 
this map). 

Extrinsic (0^) and intrinsic (X 1 ) template parame- 
ters: Extrinsic parameters are those over which p can 
be maximized efficiently, without recomputing full search 
templates for each set of extrinsic parameters under con- 
sideration (but perhaps using a small number of sub- 
templates) 30] . By contrast, maximizing p over the in- 
trinsic parameters requires computing a full template for 
each different set of intrinsic parameters in the search 
range. Searches for signals modeled by a template 
family {h(X A ) = /i(X z ,0 M )} are usually implemented 
by obtaining maxef p for each template in a discrete 
bank {h(X 1 ^, M )}, laid down only along the intrinsic- 
parameter directions. 

Minimum match MM and mismatch metric [3lj: Eqs. 



(21)-(24) of Ref. ftj. The spacing of discrete search 

banks is chosen so that at most a fraction MM is lost 

from the S/N that would be obtained with a continuous 

search bank; the corresponding loss in detection rate, for 

the same detection threshold, is a fraction MM 3 . The 

choice of the spacing is helped by considering the (full) 

mismatch metric gAB, which serves as a local quadratic 

expansion of the mismatch over all the search parameters. 

An approximation to the number of templates needed to 

achieve a given MM, computed using the metric, is given 

by Eq. (25) of Ref. [17|. 

Projected metric: Eqs. (65) and (72) of Ref. [2(j. Given 

that the search template bank has no extension along the 

extrinsic-parameter directions, it is useful to consider a 

projected metric gfj" 3 that approximates the mismatch 

already minimized over the extrinsic search parameters. 

This of, roj is still a function of the intrinsic and extrin- 
a v 

sic target parameters. The average projected metric gfj° 3 
(Eq. (75) of Ref. |2(j) is a weighted average over the ex- 
trinsic target parameters, which can be used to estimate 
the number of templates needed to achieve a given reduc- 
tion in detection rate, for a uniform distribution of tar- 
get extrinsic parameters (this reduction is proportional 

3 

to MM , where MM is the average minimum match). 

Reduced search parameter space and reduction curves 
|26l |32|: It can happen that the variety of waveforms 
spanned by an n-dimensional (search) template fam- 
ily is approximated with very high FF by an (n — k)- 
dimensional subset of the family (a reduced family) . This 
circumstance is signaled locally in the mismatch metric 
by the presence of k quasi-null directions (i.e., eigenvec- 
tors with very small eigenvalues). The integral curves 
of these directions (the reduction curves) correspond to 
sets of templates with very high match within the set, 
and map a reduced family into another. In this case, it is 
advantageous to derive the discrete search bank from a 
reduced family: ideally, one would reparametrize the full 
family using k parameters that run along the reduction 
curves, and then discard those parameters before laying 
down templates. See Ref. j2(| for a thorough discussion. 



III. SINGLE-SPIN TEMPLATE FAMILY TO 
MATCH DOUBLE-SPIN PROCESSING 
BINARIES 

This section contains the main results of this paper. 
In Sec. IIII Al we describe the PN equations for the circu- 
lar adiabatic inspiral of a double-spin binary; this tar- 
get model is used throughout this paper to represent 
physical signals. In Sec. IIII 51 we describe our proposed 
single-spin search template family, and we evaluate its 
cffectualness (which is excellent) in approximating the 
target waveforms. In Sec. IIII Gl we compare single-spin 
signals obtained at different PN orders, to argue that 
the circular adiabatic model of inspirals used in this pa- 
per gives robust predictions for the actual physical wave- 
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forms. Last, in Sec. IIIIDl we study the precessional dy- 
namics of double-spin binaries to understand which of 
its features can be represented accurately by single-spin 
systems, and which cannot. 



spin precessing binaries. The first derivative of the or- 
bital (an g ular) freq uency, up to 3.5PN order reads 

H H 13 EH El El 13 111 = 



A. Target model: double-spin precessing binaries 

Post-Newtonian calculations provide the following set 
of equations describing the adiabatic evolution of double- 
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Here, mi and m 2 are the masses of the two bodies, 
with mi > m 2 ; M = mi + m 2 is the total mass, and 
i] = mxm^/M is the symmetric mass ratio; Ljv = /ixxv 
(with ji — mim 2 /M) is the Newtonian angular momen- 
tum (with x and v the two-body center-of-mass radial 
separation and relative velocity), and Ljy = Ljv/|Ljv|; 
Si = xi m i Si and S 2 = \i m \ S2 are the spins of the 
two bodies (with Si j2 unit vectors, and < X1.2 < 1 
for BHs); je = 0.577... is Euler's constant; last, 8 = 
6 + 1987/1320+ 7cl> s /11 (with 8 an unknown arbitrary 



parameter that enters the GW flux at 3PN order [l3j . 
and lu s = 

El El El)- (Note for v3 of this paper on 
gr-qc: Eqs. (5) and (7) are now revised as per Ref. 
the parameter 8 has been determined to be 1039/4620 

00 

In Eq. Q we did not include the quadrupole-monopolc 
terms [341 . because we have already shown (see Sec. Ill 
E of Ref. |2g) that those terms do not significantly affect 
matches once these are maximized on binary parameters. 
The precession equations for the spins read |10 .18] 
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and the precession equation for Ljy is [TollTsI 



, m 2 

mi 



mi 
> ^_ 

m 2 



x L 



iV 



3 a; 1 / 3 r 
r) M 5 / 3 



(S 2 • Ljy) Si + (Si ■ Lat) S 2 



x L 



(10) 



We stop the adiabatic evolution when the binary reaches 
the minimum energy circular orbit (MECO) defined by 
Eqs. (H)-(13) of Ref. 25]. The leading-order mass- 
quadrupole gravitational waveform can be obtained from 
Eqs. (65)-(78) of Ref. J25|. Those equations for the grav- 
itational waveform, together with Eqs. (fTf> — (| 1 0|) . define 
our target model for precessing double-spin binaries. 

Using the language of Ref. 25| , precessing binaries of 
spinning BHs are described by the four basic (intrinsic) 
parameters mi, m 2 , Si = |Si|, S2 = | S2 1 , by three local 
(intrinsic) parameters describing the relative orientation 
of the spins with respect to the angular momentum at a 
fiducial frequency (see Table I and Fig. 4 of Ref. [25|), 
and by five directional (extrinsic) parameters describing 
the relative orientation of the binary and the detector 
(see Table I of Ref. 25]). The waveforms depend also (if 
trivially) on two other extrinsic parameters, the initial 
phase $0 and the time of arrival to; depending on the 
context, we shall at times omit these when counting the 
number of extrinsic parameters. 



B. Search template family: single-spin binaries 

As discussed in SeaQ] the results of previous investiga- 
tions 0, |35l |2jJ |2(| suggest that the gravitational wave- 
forms emitted by double-spin binaries with comparable 
component masses can be approximated (at least for the 
purpose of detection with ground-based interferometers) 
by waveforms computed by neglecting spin-spin effects, 
and by assigning the total spin of the binary to a single 
BH. Thus, in this paper we examine the detection per- 
formance of the single-spin search family obtained from 
Eqs. |(T||- H1U|) by setting S 2 = 0. The simplified equations 
are 
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where M s = mi, + m 2s , 77,, = mi s TO 2s /Mj and Si s = 
Xis mf s Sis- The "s" subscript stands for search param- 
eters. In Eq. (| 1 1 p . "PN corrections" denotes the terms 
in Eqs. _ O that do not depend on the spins. The 
leading-order mass-quadrupole gravitational waveform is 
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FIG. 1: Distribution of FFs for lower-mass (M < 15M ) 
binary configurations. See the caption to Table |T| for an ex- 
planation of how the FF distributions were obtained. 



given by Eqs. (65)-(78) and (11)-(13) of Ref. 25], after 
setting the spin of the lighter body to zero. This com- 
pletes the definition of our single-spin search template 
family, which is parametrized by the four intrinsic pa- 
rameters M s , rj s , xis, and Ku = Si s - Lat, and by five ex- 
trinsic parameters that describe the relative orientation 
of the detector and the binary (see Sec. Ill C of Ref. [2ti|p. 
The maximization of the overlap with respect to the ex- 
trinsic parameters can be performed semi- algebraically, 
in two steps, as described in Sec. IV of Ref. |26| . 

We note that the simplified Eqs. (| 1 1 11 (| 1 3|) are exactly 
equivalent to the full Eqs. l(T|)- (|10[) in the two limits 
mentioned in Sec. I: for equal masses, if spin-spin ef- 
fects are neglected; and for mi 3> m 2 (i.e., small rj). 
To test the effectualness of the single-spin search tem- 
plates in matching our target signals for binary config- 
urations with both comparable and unequal masses, we 
computed the FF (i.e., the match, maximized over the 
intrinsic and extrinsic search parameters; see Sec. II) for 
target binaries with two maximal spins, and with masses 
(mi+ma) = (3+3)M , (6+3)M , (6+6)M , (9+3)M , 
(12 + 3)M , (1O+1O)M , (15 + 1O)M , and (15 + 15)M . 
In the (12 + 3)M case we also considered a target binary 
with xi = 1/16 and % 2 = 1, i.e. the two objects possess 
equal magnitude of spins, while the less massive one is 
maximally spinning. Search and target signals were al- 
ways obtained at 2PN order, and the computation of the 
FF was repeated for 100 or 500 (for the lighter binaries) 
randomly generated configurations of the target-signal 
local parameters L^r, Si and S 2 , assuming uniform and 
independent angular distributions. The directional pa- 
rameters of the target signals were fixed to arbitrary val- 
ues without loss of generality, since for the purpose of 



6 





(3 + 3)M Q 


(6 + 3)M 


(6 + 6)M 


(9 + 3)M 


(12 + 3)M 


templates: 


with spin 


nospin 


with spin 


with spin 


no spin 


with spin 


with spin 


equal-spin target 


FF > 0.99 


95% 


31% 


74% 


98% 


59% 


90% 


95% 


84% 


FF < 0.99 


5% 


69% 


26% 


2% 


41% 


10% 


5% 


16% 


FF < 0.95 


0% 


38% 


3% 


0% 


25% 


1% 


0% 


0% 


lowest FF 


0.9085 


0.7042 


0.9119 


0.7250 


0.6391 


0.8945 


0.9734 


0.9684 


FF 


> 0.989 


> 0.938 


> 0.986 


> 0.987 


> 0.934 


> 0.989 


> 0.990 


> 0.990 





(10 + 1O)M 


(15 + 10)M© 


(15 + 15)M 


templates: 


with spin 


nospin 


with spin 


no spin 


with spin 


no spin 


FF > 0.99 


100% 


29% 


98% 


22% 


100% 


30% 


FF < 0.99 


0% 


71% 


2% 


78% 


0% 


70% 


FF < 0.95 


0% 


34% 


0% 


46% 


0% 


31% 


lowest FF 


0.9754 


0.7142 


0.9691 


0.7138 


> 0.99 


0.7546 


FF 


> 0.990 


> 0.945 


> 0.990 


> 0.936 


> 0.990 


> 0.957 



TABLE I: Summary of FFs between the single-spin search template family and the double-spin target model. The numerical 
maximization procedure is stopped whenever a FF > 0.99 is achieved. The upper table shows results for lower-mass binaries 
(M < 15M ); the lower table for higher-mass binaries (M > 2OM ). In the first three rows of each table, we list the percentage 
of systems yielding FFs > 0.99, < 0.99 and < 0.95, in a population of 100 target systems [500 for (mi + 7712) = (6 + 3)M 
and (9 + 3)A/ binaries] with maximal spins and random, uniform distributions of initial spin and detector orientations (local 
parameters). In the fourth row we list the lowest FFs found among the population; in the last row, we list the average FFs 
[when a FF > 0.99, we use 0.99 in computing the average.] The distribution of the FFs for selected mass configurations is also 
histogrammed in Fig. Q The target and search waveforms are computed by starting the integration of the equations of motion 
at an instantaneous GW frequency of is 60 Hz and 40 Hz for upper and lower tables, respectively. For some mass configurations 
we show also the FFs for nonspinning templates (i.e., single-spin templates where xis was set to zero), and for (12 + 3)M 
binaries, we show the FFs for a target configuration with xi = 1/16 and X2 = 1 (i.e., the 5*2 is maximal and Si = S%). 



computing FFs they are degenerate with respect to the 
local parameters (see Sec. VI A of Ref. [HI). 

The results of our tests are listed in TableHl and plotted 
in Fig. ^ For comparison, in Table [I] we include also 
some FFs computed for the nonspinning search templates 
obtained by setting \is to zero in Eqs. (fTT|> — (|T5f) . Our 
numbers support our conjecture about the effectualness 
of the single-spin search family. More specifically: 

• Spin-spin effects are not important for higher-mass 
binaries such as (15 + 15)Af , (15 + 10)Af Q and 
(10 + 10)Af Q , where FFs are consistently very high; 
however, spin-orbit effects cannot be neglected, as 
shown by the relatively low FFs achieved by the 
nonspinning search family. 

• For lower-mass binaries such as (3 + 3)Af© and 
(6 + 6) Mo, FFs are also very high, with few ex- 
ceptions: thus, although in these binaries spin-spin 
effects can accumulate over many GW cycles within 
the band of good detector sensitivity, they rarely 
become comparable to spin-orbit effects. 

• For low-mass-ratio binaries such as (12 4- 3)M© 
(r] = 0.16), FFs are high, since the spin of the heav- 
ier object dominates the precessional dynamics. If 
we reduce the magnitude of Si so that Si = S2, 
the resulting FFs become lower, because the dy- 
namics deviates farther from both the single-spin 



and equal-mass limits. 

• The worst FFs are obtained for (6 + 3)M© and (9 + 
3)Mq binaries, which have rather low total masses, 
and intermediate mass ratios (thus, they sit halfway 
between the two single-spin equivalence limits). In 
this case, double-spin effects cannot be reproduced 
with accuracy by single-spin systems (in Sec. Ill D 
we shall examine in more detail what is happening 
there). Note however that this happens only for 
a limited number of angular configurations, so the 
average of the FF over the sampling is still very 
high. 

The range of search-template parameters needed to 
yield the high FFs discussed above extends beyond values 
that would be physical for a real single-spin binary, with 
■q s > 0.25 and \is > 1- This is to be expected: consider, 
for instance, that in the equal-mass case the equivalence 
between the simplified and the full equations implies val- 
ues of xi s up to 2. In Fig. [2] we show the parameters of 
the best-fit search templates corresponding to target sig- 
nals with the test masses examined above [augmented by 
(15 + 3)M Q , (12 + 6)Af , and (15 + 12)Af©]. As shown in 
the top panel, the search-template images of target sig- 
nals with the same masses but different local parameters 
are spread around the nominal (M s , r] s ) values (indicated 
by the end of the thin lines, and always enclosed within 
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FIG. 2: Location in the (intrinsic) search parameter space (M s , rj a> xist «is) of the best-fit templates for target signals with 
(mi+mj) = (3 + 3)M , (6 + 3)M Q , (9 + 3)Af , (12 + 3)M , (15 + 3)M Q , (6 + 6)Af , (12 + 6)Af , (1O + 1O)M , (15 + 1O)M , 
(15 + 12)Me, an d (15 + 15)M©, with maximal spins, and with random angular distributions of the initial Ljv, Si, S2. Dots are 
denser for the (6 + 3)M and (9 + 3)M configurations, for which more FF were computed. In the (M s , rj s ) scatter plot (on top), 
the dashed contour encloses the region obtained by setting M a = M and rj s — rj, and by taking (mi , m^) € [3, 15]M x [3, 15]M . 
The dotted and dashed line, drawn somewhat arbitrarily, encloses a possible template bank boundary, used in Sec. II Vl to estimate 
the number of templates necessary to search for double-spin binaries in this mass range. The labels identify the search template 
clusters corresponding to each target mass configuration, and they are connected to the nominal projection point obtained by 
setting M a = M and rj s — rj. 



the dashed contour). The uncertainty in target parame- 
ter estimation induced by this spreading is discussed in 
Sec. [V] In the same panel, the dotted-dashed line en- 
closes the template-bank boundary used in Sec. IIVI to 
estimate the number of templates necessary for a search 
of systems with masses (mi,m 2 ) € [3, 15]M© x [3, 15]M Q . 
The bottom panels show the range achieved by the search 
parameters xis and ki s , which is comparable to the 
range of the analogous target parameters, |Stot|/wf and 

K-tot = Stot ' Lat. 



C. On the robustness of waveforms across PN 
orders 

In the previous section we have established that single- 
spin waveforms are good approximations for double-spin 



waveforms, at least within the mass range under con- 
sideration. However, whether double-spin waveforms are 
representative of actual physical signals is an entirely dif- 
ferent question, which hinges on the validity of the cir- 
cular adiabatic approximation, but also on the robust- 
ness of the waveforms under change of PN order: if the 
waveforms change substantially with increasing order, we 
should suspect that the description of the physics is in- 
complete without higher-order terms yet to be computed. 

Studying the robustness of double-spin binaries is tech- 
nically difficult, since it means computing the FF be- 
tween two template families (of different PN order) with 
seven intrinsic parameters. This entails the delicate 
numerical maximization of a seven-parameter function 
whose evaluation is relatively costly. Instead, we choose 
to perform our study on single-spin waveforms, and then 
argue that the results should transfer to double-spin 
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(N + k,N) 



(STiv+fe, STat) for (10 + 10) M 
ki = 0.9 ki = 0.5 



binary, M = 2OM , r\ = 0.250 
Ki = -0.5 



ki = -0.: 



(1,0) 

(1.5,1) 

(2,1.5) 

(2.5,2) 

(3,2) 

(3,2.5) 

(3.5,3) 



0.3136 (0.6688) [16.1,0.25,0.00, 0.00] 

0.3123 (0.5922) [23.7,0.25,0.00, 0.00] 

0.7124 (0.9823) [19.2,0.25,1.00, 0.99] 

0.2851 (0.8702) [18.8,0.25,1.00, 0.91] 

0.9604 (>0.99) [20.2,0.25,0.88, 0.79] 

0.2848 (0.9846) [18,8,0.25,1.00,-0.31] 

0.9316 (>0.99) [20.2,0.25,1.00, 0.90] 



0.3136 
0.2676 
0.7222 
0.3099 
0.9743 
0.2898 
0.9475 



(0.6688) 
(0.5137) 
(0.9877) 
(0.9206) 
(>0.99) 
(>0.99) 
(>0.99) 



0.3136 (0.6688) 

0.2306 (0.4860) 

0.8545 (0.9886) 

0.4166 (>0.99) 

0.9915 (>0.99) 

0.4027 (0.9823) 

0.9749 (>0.99) 



[16.1,0.25,0.00,-0.00] 
[29.2,0.25,0.00,-0.00] 
[19.3,0.25,1.00,-0.66] 
[20.4,0.25,1.00, 0.81] 
[19.7,0.24,0.88,-0.23] 
[20.4,0.25,1.00,-0.99] 
[19.7,0.25,1.00,-0.99] 



0.3136 (0.6688) 

0.2160 (0.4543) 

0.8601 (>0.99) 

0.4682 (>0.99) 

0.9805 (>0.99) 

0.4634 (0.9740) 

0.9744 (>0.99) 



(N + k,N) 



(STjv+fc, STat) for (15 + 1O)M binary, M = 25M , 17 = 0.247 
0.9 «i = 0.5 Ki = -0.5 



Kl 



-0.9 



(1,0) 

(1.5,1) 

(2,1.5) 

(2.5,2) 

(3,2) 

(3,2.5) 

(3.5,3) 



0.3124 (0.6030) 

0.2784 (0.4994) 

0.5810 (>0.99) 

0.2558 (0.8525) 

0.9106 (>0.99) 

0.2520 (0.9148) 

0.9264 (>0.99) 



[19.1,0.25,0.00,0.00] 
[31.0,0.25,0.00,0.00] 
[24.8,0.23,1.00,0.98] 
[22.9,0.23,1.00,0.95] 
[24.8,0.24,0.88,0.73] 
[27.0,0.25,1.00,0.85] 
[24.9,0.24,1.00,0.90] 



0.3124 (0.6030) 

0.2732 (0.4684) 

0.8038 (>0.99) 

0.3296 (0.9280) 

0.9392 (>0.99) 

0.2942 (>0.99) 

0.9528 (>0.99) 



0.3124 (0.6030) [19.1,0.25,0.00,-0.00] 

0.2466 (0.3913) [41.9,0.25,0.00,-0.00] 

0.8644 (>0.99) [25.4,0.24,1.00,-0.59] 

0.5022 (>0.99) [25.5,0.24,1.00, 0.80] 

0.9858 (>0.99) [25.1,0.23,0.89,-0.21] 

0.4552 (>0.99) [25.7,0.24,1.00,-0.59] 

0.9769 (>0.99) [25.3,0.24,1.00,-0.36] 



0.3124 (0.6030) 
0.1896 (0.3491) 
(>0.99) 
(>0.99) 
(>0.99) 
(>0.99) 
(>0.99) 



0.9067 
0.5921 
0.9650 
0.5333 
0.9839 



(TV + k,N) 



(ST N+k , STat) for (15 + 15)M 
0.9 «i=0.5 



binary, M = 3OM , t] = 0.250 
m = -0.5 



-0.9 



(1,0) 

(1.5,1) 

(2,1.5) 

(2.5,2) 

(3,2) 

(3,2.5) 

(3.5,3) 



0.2710 (0.5158) [22.5,0.25,0.00, 0.00] 

0.2694 (0.4050) [38.3,0.25,0.00, 0.00] 

0.7619 (>0.99) [31.1,0.25,1.00, 0.89] 

0.3403 (0.9086) [28.0,0.24,0.92, 0.99] 

0.9330 (>0.99) [28.8,0.24,0.90, 0.82] 

0.2926 (>0.99) [28.1,0.25,1.00,-0.58] 

0.9265 (>0.99) [30.4,0.25,0.90, 0.71] 



0.2710 
0.2145 
0.8613 
0.3856 
0.9360 
0.3549 
0.9456 



(0.5158) 
(0.3644) 
(>0.99) 
(0.9237) 
(>0.99) 
(>0.99) 
(>0.99) 



0.2710 (0.5158) [22.5,0.25,0.00,-0.00] 
0.2435 (0.3155) [49.6,0.25,0.00,-0.00] 
0.9018 (>0.99) [30.6,0.25,1.00,-0.43] 
0.5372 (>0.99) [30.4,0.25,1.00, 0.95] 
0.9641 (>0.99) [30.3,0.24,0.88, 0.17] 
0.4893 (>0.99) [31.0,0.25,1.00,-0.67] 
0.9814 (>0.99) [29.1,0.25,1.00,-0.30] 



0.2710 (0.5158) 
0.1855 (0.2797) 
(>0.99) 
(>0.99) 
(>0.99) 
0.5498 (>0.99) 
0.9831 (>0.99) 



0.8946 
0.6216 
0.9612 



TABLE II: Test of robustness of the PN adiabatic waveforms STat (defined in Sec. 1111 Al across PN orders, for (mi + 1712) = 
(1O+1O)M , (15 + 10)M Q and (15 + 15)M . We set X i = 1, X2 = 0, and «i = -0.9, -0.5,0.5, and 0.9. The matches quoted at 
the beginning of each column are maximized only with respect to the extrinsic parameters to and <3>o- In parentheses, "(...)", we 
give the matches maximized over all the parameters of the lower-order family (i.e., the fitting factors FF for the target family 
STAf+fc as matched by the search family STat, evaluated at the STAr+fc intrinsic parameters indicated). In brackets, "[...]", we 
give the parameters M, r/, xi, and ki (or M and r\ at 1PN and 1.5PN orders when spin terms are absent) at which the FF 
is achieved. The detector is set perpendicular to the initial orbital plane, and at 3PN and 3.5PN order we set 9 = 0; in all 
cases the integration of the equations of motion starts at an instantaneous GW frequency of 40 Hz. [See Refs. [17, 25] for a 
discussion of why for some mass combinations the 2.5PN model differs so much from the other orders.] 



waveforms because the search and target families are 
close. Thus, our study is very similar to the Cauchy con- 
vergence test of Ref. j26(, except for the choice of masses: 
here we focus on binaries at the higher-mass end of our 
range, since these systems are expected to have stronger 
higher-order PN effects within the frequency band of 
good interferometer sensitivity. 

For (mi + m 2 ) = (10 + 1O)M , (15 + 10)M Q , and 
(15+15)A/ , we list in TablelTTlthe matches across PN or- 
ders, maximized only on to and $0; the numbers in paren- 
theses, "(■•■)", give the FF for the higher-order family as 
matched by the lower-order family, and the numbers in 
brackets, "[■■■]", give the intrinsic parameters where the 
FF is attained. The tests are performed for xu = 1, for 



different values of «i s , and for a GW detector in a direc- 
tion (with respect to the binary) perpendicular to normal 
vector of the initial orbital plane, which should be rep- 
resentative of the generic effects of precession. The high 
FFs obtained between the 2PN and the higher-order fam- 
ilies suggest that the 2PN model is already representative 
of the variety of waveforms expected from actual sources; 
on the other hand, the lower direct matches (and the bi- 
ased values of search parameters at the FF) suggest that 
the family of the highest available order should be used 
for source parameter estimation. It would be worthwhile 
to evaluate the FF between the double-spin (and indeed, 
single-spin) adiabatic model, and nonadiabatic models 
based on resummed PN equations 0, 0, |3(J , espe- 
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FIG. 3: Relative change of the opening angles as function of 
8ls for a (6 + 3)Af© binary, with xtot = 0.4, |Lj = |Ljv| = 
r?M 5/3 u~ 1/3 , and ui = 2ix x 30 Hz. The change shown cor- 
responds to a 10% increase in |Stot|- The solid and dashed 
curves refer to 9l and 9s, respectively. 



cially when these predict the end of the inspirals within 
the band of good interferometer sensitivity. 



D. Some features of the dynamics of double-spin 
binaries 

In this section we study the precessional dynamics of 
double-spin binaries, with the purpose of building a phys- 
ical understanding of the matching performance of single- 
spin templates observed in Sec. IIII Bl in particular, we 
wish to identify what features of double-spin dynamics, 
absent in single-spin systems, lead to the low FFs seen 
for lower-mass binary configurations. 

From Eqs. (|8|)- i|Ttj|) . we see that the preces- 
sion of double-spin binaries preserves both J and 
|L| at timescalcs shorter than the radiation- reaction 
timescale — at which |L| decreases steadily. Even at 
longer timescales, as recognized by ACST for single-spin 
binaries and further tested by BCV2 for generic double- 
spin binaries, for the vast majority of configurations the 



direction of the total angular momentum remains almost 

constant (J ~ 0); this behavior is known as simple pre- 
cession. 

For single- spin binaries, or for equal-mass binaries if we 
ignore the spin-spin interaction, the angle between L and 
Stot (@ls) remains fixed all through evolution [according 
to Eqs. (JSJ) (|10fl ]: f° r simple precession, this implies that 
the angle between L and J (9l) must increase, and that 
the angle between S to t and J (9s) must decrease — both 
do so monotonically, at the radiation-reaction timescale. 
In summary, in these binaries the orbital plane precesses 
while its inclination increases slowly and monotonically. 

In Ref. [2lJ , Apostolatos investigated the effect of spin- 
spin coupling on the dynamical evolution of equal-mass, 
equal-spin BH-BH binaries. He obtained analytical solu- 
tions for the opening-angle products Si. 2 ■ J and Si • S2 
to first order in S/L, where J is the total angular mo- 
mentum, S is Si or i§2) an d L is the orbital angular mo- 
mentum (if mi ~ m2, then throughout all the inspiral 
S -C L [13 ). The main feature identified by Aposto- 
latos was that the orbital plane not only becomes slowly 
more inclined (at the radiation-reaction timescale), but 
the spin-spin interaction also causes a nutation; namely, 
an oscillation of the orbital inclination angle (9^) at the 
timescale of the spin-spin interaction. 

In the following, we shall relax the assumption that the 
BH masses are equal, and we shall investigate what the 
consequences are on the evolution of Si^ • J and Si ■ S2. 
To simplify our notation we fix mi + 777,2 = 1, and we 
introduce the parameter 6 = mi — 777,2, which describes 
the deviation from the equal-mass case. We keep only 
terms up to linear order in S. We then have (assuming 
as always m\> 1112) 



1 



1-5 



mi = 



m 2 



S 1 =m\ Xl ~-(\ + 25)xi 



S2 = m 2 X2 - 7 (1 - 2<5) X2 



(14) 
(15) 
(16) 



Inserting these definitions into Eqs. JSJ and (JSJ leads to 



J 



dt 



w2 s 
T 2 



(S 2 -J) = 



spin— orbit 

(7-65) - ( 1 + :-! 

spin— orbit 



j • s ; 



L 2 

spin— spin 




J • (Si x S 2 ) 



J - (S 2 x Si) 



J -(Six S 2 ) 



(17) 



(18) 



(19) 
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Following Apostolatos, in deriving these equations we 
have assumed that S <C L, L = J [1 + 0(S/L)], and that 
the direction of the total angular momentum remains al- 
most constant during evolution. In Eqs. (|17[) ~ (|19[) we 
have separated the terms due to spin-orbit and spin-spin 
interactions. 

According to Eq. JT^J, in the equal-mass case (5 = 0) 
without spin-spin effects, the angle between Si and S2 is 
constant; generically, however, the spin-spin and even the 
spin-orbit interactions can cause that angle (and hence 
the magnitude of Stot) to change. As first observed by 
ACST and further investigated by Apostolatos |21|. this 
variation in |Stot| drives the nutation of the orbital plane: 
oscillations are superimposed to the monotonic evolution 
of the angles between L and J (9l) and between S to t 
and J (9 s) |2l|, as can be understood from the following 
simple argument. Recall that on timescales shorter than 
the radiation-reaction timescale, |L| and J are conserved; 
using 8 to denote the change in the dynamical variables 
incurred during such a time, we write 



S\J\ 2 =0, 8\L\ 2 = 



to get 



2<J(L-S tot ) = -£(|S tot | 2 ). 
Using 83 = we then have 

5\L = -^-M\S tot \ 2 ) for X L LJ 



2LJ 



LJ 



(20) 



(21) 



(22) 



and 



<5A5 = -fe|^(|Stot| 2 ) for A s = %^. (23) 



Stot J 



Thus, when | Stot | oscillates, the opening angles Oh and 
9s oscillate as well; in fact, we have 



SOt 



1 



2|Lx S to t 1 
L • S t ot 



2|Stot| 2 |Lx S to t 



<5(|Stot| 2 ), (24) 
<5(|Stot| 2 ), (25) 



which suggests that, for the same variation in |Stot| 2 , 
the nutation is most significant when L and S are ei- 
ther nearly aligned or antialigned. In Fig. |31 we plot 
the relative changes in 9l and 9s as functions of the an- 
gle between L and S (9ls = 9l + 9s), choosing a fixed 
positive <5(| Stot | ) - The change is always positive for 9l, 
while it can be negative for 9s, if #ls > 90°. In addition, 
the relative changes diverge near Ols ~ 0° or 180°. These 
features follow straightforwardly from Eqs. I|24|) and (|25|l . 

Spin-spin effects. When spin-spin effects are included, 
the angle between Si and S2 oscillates according to the 
second term on the right-hand side of Eq. (|19|) . How- 
ever, as evidenced by the FF results for equal-mass bi- 
naries given in Sec. Ill B, the amplitude of these oscil- 
lations does not seem to be very large, and the nutation 



of the orbital plane does not complicate significantly the 
waveforms, at least as evaluated at the leading mass- 
quadrupole order. 

Spin-orbit effects. From the first term on the right- 
hand side of Eq. 1|19|) , we see that even in the absence of 
spin-spin effects, unequal masses (i.e., 8 ^ 0) can cause 
the evolution of Si • S2, and therefore drive the nutation 
of the orbital plane. Spin-orbit effects, which come in at 
a lower PN order, can sometimes be more significant than 
spin-spin effects, especially for binaries with intermediate 
mass ratios, such as (7111+1112) = (6+3)A/q and (9+3)Af 
binaries; indeed, these effects could explain the lower FFs 
found in Sec. Ill B for those systems. 

Examples of combined spin-spin and spin-orbit effects 
in double-spin binaries leading to oscillations in |S to t|, 
and therefore 9^ and 9s, are shown in Fig. 0] for systems 
with (mi + 1112) = (6 + 3)M Q (8 = 1/3), and with ini- 
tial local parameters such that FF > 0.99 (on the left) 
and FF ~ 0.94 (on the right). The nutation behavior 
evident in these figures is well described by the approxi- 
mated equations l|24l) and l|25|l . For comparison we show 
(as continuous lines) also the evolution of the analogous 
quantities in the best-fit single-spin configurations; for 
these, the opening angles 9l and 9s evolve monotonically. 
Lower overlaps seem to correspond to initial conditions 
for which nutation is rather significant and overwhelms 
the underlying monotonic evolution. In Fig. we show 
the percentage of configurations with FF < 0.99 (light 
pattern) and FF < 0.97 (dark pattern) as a function of 
the initial Xj,s = Ljy ■ Stot (evaluated at the starting 
frequency 38]), for (6 + 3)M and (9 + 3)M© binaries. 
The plot suggests that lower FFs are more likely to occur 
when the initial L • S to t — ±1, which is consistent with 
the fact that nutation is most significant when 9ls ~ 0° 
or 180° EHI. 

Figure shows also an asymmetry in the distribution 
of low FFs, which arc denser when 9ls < 90°. Currently 
we do not have a clear understanding of this behaviour 
|40j . It is worth pointing out that this asymmetry and, 
more in general, the low FFs observed could be due to 
other features of double-spin dynamics that cannot be 
reproduced by single-spin systems, but that are difficult 
to dig out by analysis or numerical experiment, in the ab- 
sence of a full analytical solution to the precession equa- 
tions. Moreover, some of the low FFs might be due to 
shortcomings in our numerical optimization procedure in 
cases where the match surface in the search parameter 
space has an especially convoluted geometry. 



IV. TEMPLATE SPACE AND NUMBER OF 
TEMPLATES 



In this section we estimate the number of single-spin 
templates necessary to search for double-spin signals with 
single masses in the [3, 15]M Q range. To do this, we com- 
pute the average projected metric gf™ 3 in the (M s ,i] s ) 
region delimited by the dashed and dotted contour of 
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FIG. 4: Evolution of the opening angles 6l and 6s, and of the total-spin magnitude Stot (all plotted as dashed lines) for 
double-spin target systems yielding FF > 0.99 (left column) and FF ~ 0.94 (right column) when matched by single-spin 
templates; the target system has (mi + m.2) = (6 + 3)Mq. For comparison, the solid lines show the evolution of the analogous 
single-spin quantities [6l, 8s ~ arccos(Si s ■ J), and Si s ] for the best-fit single-spin systems. 



Fig. following the procedure described in Sec. VI of 
Ref. |26j. We notice the presence of reduction curves 
connecting (roughly) the segment xis G [0, 2],ki s = 
to the entire (xisjMs) plane in the search template 
space. Thus, we select a 3-D reduced template space 
corresponding to (M s ,rj s ) within the quadrilateral with 
vertexes (15M Q ,0.08), (3.25M Q , 0.275), (32.5M , 0.23), 
(24.5M Q , 0.385), to xu G [0,2], and to ku = 0. Addi- 
tional subfamilies might be needed to deal with certain 
singularities that we observe in the reduction curves as 
Ku gets close to ±1, but our selection should already 
give us an acceptable idea of the number of necessary 
templates, which is computed according to 



templates 



/ 



dM s dr] s dxis 



2y/(l -MM)/3 



(26) 



where the primed indices i', j' run through M S} rj s , and 
Xisj the metric is averaged over 1,000 sets of target ex- 
trinsic parameters. The integral is carried out by evalu- 
ating the projected metric (a computationally expensive 
operation) at 80 points within the integration region, and 



filling it by natural neighbor interpolation |4l| . The final 
result is Mompiatos ^ 320,000 for MM = 0.98 (not includ- 
ing a reduction mismatch of ~ 0.01 incurred along the 
reduction curves |2(j). Given the uncertainties implicit 
in the numerical computation of the metric, in the inter- 
polation, in the choice of the reduction curves, and in the 
actual placement of the templates in the bank, this num- 
ber should be understood only as an order-of-magnitude 
estimate. Most of the templates, by a factor of many, 
come from the lower part of the integration region (i.e., 
from the lowest ?y s for any given M s ); about 100,000 out 
of 320,000 come from the region with Xis < !• 



V. ESTIMATION OF BINARY PARAMETERS 

Since the single-spin template family contains only a 
subset (of lower dimensionality) of all possible double- 
spin waveforms, we cannot expect to obtain estimates of 
all the physical parameters of double-spin systems from 
a single-spin-template search. The most straightforward 
way to recover those parameters would be to perform a 
search using double-spin templates after single-spin tern- 
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plates have yielded a detection; such a follow-on search 
may be computationally feasible, since double-spin tem- 
plates will then be applied only to the data stretches 
that have been established to contain signals. However, 
it is still meaningful to analyze the parameter-estimation 
performance of the single-spin template family, since any 
constraints on source parameters will decrease the size of 
the double-spin template bank necessary for the follow-on 
search, lowering the computational cost even further. In 
addition, this analysis can offer further useful insight into 
the FF map of double-spin into single-spin waveforms. 

Double-spin target waveforms are parametrized by 
twelve parameters (seven intrinsic, and five extrinsic), 
while single-spin search templates are parametrized by 
nine (four intrinsic, M s , ?7 S , Xls, & n d k\ s ; and five ex- 
trinsic). Thus, maximizing the match over the search 
parameters induces a map from the 12-parameter target- 
signal space into the 9-parameter template space. The 
inverse map takes each point in the template space into a 
3-dimensional manifold in the target-signal space, whose 
size indicates the extent to which the physical parameters 
of the double-spin binary can be constrained. Unfortu- 
nately, evaluating the size of the inverse image requires 
computational resources well beyond what is currently 
available to us. Furthermore, such a procedure can only 
be meaningful after statistical error has also been taken 
into account: because the detector output contains also 
noise, the template parameters at which we obtain the 
maximal correlation between template and data (the ac- 
tual projection point) will differ, by a random statisti- 
cal error, from the parameters at which the correlation 
between template and signal is highest (the theoretical 
projection point) 0- We leave the quantitative study 
of statistical errors to a forthcoming paper • 

In this paper we take a semi-quantitative approach. 



We ignore the extrinsic parameters, so we study the map 
of seven parameters into four; then, we identify a num- 
ber of intrinsic parameters of the double-spin binary (our 
target observables) , and we explore how well we can es- 
timate their values using functions of the four intrinsic 
parameters of the best-fit single-spin template (our esti- 
mators). We use three general criteria in the choice of 
target observables and estimators: 

Consistency. The target observables and estimators 
should coincide when the double-spin system is dy- 
namically equivalent to a single-spin system (ac- 
cording to the criteria spelled out in Sec. I). 

Robustness. The definitions of the target observables 
and of the estimators should be independent of the 
detector noise curve; equivalently, the target ob- 
servables should be (almost) conserved quantities 
(we already know that the four template parame- 
ters, and hence the estimators, are conserved). 

Strong influence. The target observables, and there- 
fore the estimators, should have a strong influ- 
ence to the waveforms. Quantitatively, we can re- 
quire the mismatch metric to have large compo- 
nents along the direction of change of the target 
observables and estimators; it follows that the tar- 
get observables and estimators should remain es- 
sentially constant along eventual reduction curves. 

This criterion is important for two reasons: first, 
the 7-to-4-dimensional FF map in unlikely to pre- 
serve unessential features of the target space; sec- 
ond, even if the map preserved these features, sta- 
tistical error would inevitably spoil their estima- 
tion, because the associated mismatch-metric com- 
ponents are small. A third reason, contingent on 
our implementation of the FF search, is that we 
stop the maximization of the match whenever this 
reaches 0.99; this adds a dominant artificial fluctu- 
ation (roughly corresponding to the statistical error 
for a S/N of 100) to parameter estimation. 

Since our search template family possesses a family of 
approximate reduction curves, it will be generally possi- 
ble to estimate efficiently only three independent target 
observables, and four only for very high S/N. Our hope 
is then to find three target-parameter-estimators pairs 
that satisfy all three criteria, and four pairs that satisfy 
the consistency and robustness criteria. 

It is straightforward to see that using M s to estimate 
M and rj s to estimate rj (or in shorthand, M s — > M and 
•q s — > 77) satisfies the consistency and robustness crite- 
ria. We shall pay a special attention to the estimation of 
chirp mass according to M. s = M s r)s — * M. = Mrj 3 ^ 5 , 
which satisfies all three criteria (in particular, chirp mass 
is conserved very accurately along reduction curves). In 
upper part of Table IIHI we characterize the distribution 
of estimation error for M, 77, and A4, for the same sys- 
tems (with fixed masses and random local parameters) 
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(6 + 3) M Q 
(500 points) 


(9 + 3) M Q 
(500 points) 


(12 + 3)M 
(100 points) 


(10 + 10) M 
(500 points) 


(15 + 10) M 
(100 points) 


(15 + 15) M Q 
(100 points) 


M 


{M s -M)/M 
AMs/M 
l-a/3-a percentage 


+0.0232 
0.0817 
75.6%/98.0% 


+0.0021 
0.0755 
84.0%/97.4% 


-0.0066 
0.0631 
79.0%/97.0% 


-0.0351 
0.0437 
76.4%/98.6% 


-0.0206 
0.0603 
70.0%/99.0% 


-0.0579 
0.0639 
72.0%/100.0% 


V 


Vs-V 
At] s 

l-a/3-a percentage 


-0.0057 
0.0268 
72.6%/99.0% 


+0.0022 
0.0254 
84.4%/98.8% 


+0.0038 
0.0178 
81.0%/98.0% 


+0.0191 
0.0219 
76.6%/98.6% 


+0.0122 
0.0283 
74.0%/99.0% 


+0.0362 
0.0387 
74.0%/100.0% 


M 


(Ms - M)/M 
AMs/M 
l-a/3-a percentage 


-0.0004 
0.0074 
71.6%/99.0% 


+0.0015 
0.0106 
69.6%/99.2% 


+0.0021 
0.0104 
66.0%/100.0% 


+0.0055 
0.0092 
71.2%/99.2% 


+0.0033 
0.0145 
70.0%/98.0% 


+0.0142 
0.0192 
80.0%/99.0% 




i LjV'^off 1 JV ' ® c If 
V M ' 2 )s M ' 2 

A 1 / I'.V Snff \ twS of f] 
M 2 )s M ' 2 1 

l-a/3-cr percentage 


+0.0236 

0.1290 
79.4%/98.0% 


+0.0210 

0.1130 

73.4%/98.8% 


+0.0129 

0.1010 
68.0%/100.0% 


-0.0071 

0.0747 
75.8%/99.0% 


-0.0015 

0.0886 
66.0%/100.0% 


+0.0004 

0.0917 
73.0%/98.0% 


& cft 

M 


M \ 

1-a/S-a percentage 


+0.0253 

0.1920 
75.0%/98.6% 


-0.0264 

0.1480 
72.6%/98.6% 


-0.0183 

0.1050 
67.0%/99.0% 


+0.0375 

0.1170 

75.0%/98.8% 


+0.0028 

0.1860 
74.0%/100.0% 


+0.0387 

0.1460 
79.0%/99.0% 


Xtot 


(Xtot)s Xtot 

A [(xtot) s - Xtot] 
l-a/3-a percentage 


+0.0627 
0.2040 
81.6%/98.2% 


-0.0053 
0.1610 
81.8%/97.8% 


-0.0140 
0.1190 
85.0%/97.0% 


+0.0174 
0.0901 
78.8%/98.0% 


+0.0036 
0.1450 
71.0%/99.0% 


+0.0149 
0.1220 
80.0%/98.0% 



TABLE III: Systematic biases, rms deviations, and percentage of samples within ±1 and 3 deviations of the average, for six 
target-observable-estimator pairs. The mass configurations are those studied in Sec. 1111 61 



used in Sec. IIII Bl to compute FFs. Each section shows 
the estimation bias (defined as the average of the er- 
ror, and measuring a systematic displacement between 
observables and estimators that can in principle be re- 
moved), its rms deviation (measuring the intrinsic un- 
certainty in the estimation), and the percentage of esti- 
mators enclosed within 1-dcviation and 3-deviation inter- 
vals (measuring the normality of the distribution). The 
chirp mass, which satisfies all three criteria, is indeed es- 
timated with higher relative accuracy than both M and 
7]. Even after statistical errors are taken into account, 
this accuracy should be retained for M. better than for 
M and rj. The distributions of M s , r/ s , and M. s are also 
histogrammed in Fig. EI for + m 2) = (6 + 3)M and 
(10 + 10)M Q binaries. 

It is very hard to identify additional target-parameter- 
cstimator pairs that satisfy all three criteria, mainly be- 
cause double-spin binaries lack conserved quantities that 
clearly dominate the waveforms; so the two additional 
target observables that can be estimated efficiently may 
not have simple physical meanings. When spin-spin ef- 
fects are negligible, the only truly conserved quantity 
that could be interesting for our purposes is S c ff • Ljvj 



with 

[The magnitudes of the individual spins do not satisfy the 
consistency criterion, since single-spin binaries (e.g., with 
mi spinning) require |Si s | — > |Si|, while double-spin, 
equal-mass binaries require |Si a | — > |Si + Saj.] There- 
fore, we choose the target observable S e g -ILn/M 2 , which 
is conserved and closely related to the opening angle be- 
tween S c ff and Ljv , and hence to the depth of the mod- 
ulation caused by orbital precession. The estimator is 
naturally 

S c ff -Lat _ / 3m 2s \ m\ s 

J s 

where (mp), = Re [(1 ± y/1 - /2] M s (taking the 
real part becomes necessary when rj s > 0.25). As we 
see from Table IIIII this observable can be estimated 
with bias within [—0.01, +0.02], and with rms devia- 
tion 0.07-0.13. Although conserved, this observable is 
not quite constant along reduction curves. In the left 
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FIG. 6: Distribution of errors for the target observables M, r), and M, as estimated by M s , rj a and M a , for 500 double-spin 
binaries with (mi + m-z) = (6 + 3)Mq and (10 + 1O)M0, maximal spins, and uniform distributions of local parameters. The 
Ms distribution has the smallest bias and dispersion; the M s and r) s distributions have much larger dispersion and are skewed 
in opposite directions (as needed to reduce the dispersion of M s = MaTj^ 5 ). 



panels of Fig. |7| we plot the distribution of the pairs 
(S eS -t N /M 2 , [S cff -t N /M%) for ( mi +m 2 ) = (6 + 3)M 
and (10 + 10)M© binaries. In each case we conclude 
that the target observable and the estimator are strongly 
correlated; however, the dispersion is noticeably smaller 
for (10 + 1Q)Mq than for (6 + 3)M binaries; this must 
be because for equal-mass binaries only spin-spin effects 
can cause differences between double-spin and single-spin 
waveforms. 

In the light of Eq. (|10f) , we choose the third target ob- 
servable as | S c fi | /M, which measures the instantaneous 
angular precession frequency divided by lu 2 . This quan- 
tity is conserved only for single-spin binaries and for 
equal-mass binaries with negligible spin-spin effects, so it 
does not completely satisfy the robustness criterion. In 
our study, we use the value of |S e ff|/M at the initial fre- 
quency (40 or 60 Hz) from which the equations of motion 
are integrated. [However, it would be more reasonable to 
evaluate it at the frequency at which the detector is most 
sensitive, or to weight its values at different frequencies 
according to the detector noise spectrum.] The estimator 



is 

To make the observable and the estimator dimensionless, 
we divide both by M ; TablelTTTIthen shows that |S cff \/M 2 
can be estimated with bias within [—0.03, +0.04], and 
rms deviation 0.11-0.19. Despite the apparent physi- 
cal meaning of [|S e ff|/M] s , this observable is also not 
conserved well along reduction curves. Thus, we might 
as well use a more familiar target observable, Xtot = 
S tot /M 2 , as estimated by [xtot] s = Xism\ s /M 2 . This 
pair satisfies the consistency criterion, and it changes 
through the inspiral at a level similar to S c ff/M 2 . Ta- 
ble shows that xtot can be estimated with bias within 
[-0.005, +0.06], and rms deviation 0.10-0.20. In the cen- 
ter and right panels of Fig. we plot (target, estimator) 
distributions for these two spin observables, again for for 
(mi + m 2 ) = (6 + 3)M Q and (10 + 1O)M binaries. It 
is clear from the plots that the accuracy of estimation is 
poorer than for S c ff ■ Lat/M 2 — > [S c ff ■ Ljv/M 2 ] s ; again, 
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the dispersion is smaller for the (10+10)M Q binaries. For 
(6 + 3)M binaries, the worse accuracy can be attributed 
in part to the fact that these two target observables are 
not as well conserved as S c g • /M 2 during evolution. 

VI. CONCLUSIONS 

As originally pointed out by ACST , the dynamics 
of double-spin precessing binaries become equivalent to 
the dynamics of single-spin binaries (at least for the pur- 
pose of computing gravitational waveforms at the lead- 
ing mass-quadrupole order) in two limits: equal masses, 
when spin-spin effects can be neglected (then Si — > Stot), 
and very different masses (then Si tends to the spin of 
the heavier body). Building on this observation, on the 
results of Refs. 0, 0, EBf . and on the (justified) as- 
sumption that spin-spin effects contribute mildly to the 
PN binding energy and GW flux of the binary for mass 
configurations of interest to ground-based GW interfer- 
ometers, we conjectured that single-spin templates (as 
defined in Sec. lIII B|l can be used effectually to search for 
double-spin precessing binaries with such masses. 

We tested our conjecture by evaluating the FF between 
the single-spin and double-spin families, and we found 
confirmation in the very high FF values [see Tabled and 
Fig. ^ for equal-mass binaries of both low and high total 
masses. FFs were high also for unequal-mass binaries, 
except for few initial spin configurations. As discussed 
in Sec. HID, for those configurations the evolution of the 
opening angles between J and Ljv and between J and Si,2 
seem to contain large oscillations, induced by spin-spin 
and non-equal-mass effects, that cannot be reproduced 
sufficiently well by single-spin systems. 

The region in the single-spin parameter space 
needed to match double-spin binaries with (011,7713) = 
[3, 15]M x [3,15]M Q is shown in Fig. [3 Using the 
LIGO-I design sensitivity, we counted (very roughly) as 
~ 320,000 the number of templates required to yield a 
minimum match of 0.97. The number of BCV2 templates 
needed for a similar mass range is somewhat larger. More 
generally, with respect to the detection template families 
introduced in Refs. |U lU IH |H ^ , the advantage of 
the quasi-physical single-spin family suggested in this pa- 
per is the possibility of estimating the parameters of the 
source. In Sec.|V]we computed the systematic errors that 
would affect the measurement: the total mass M could 
be estimated with a fractional bias within [—6%, +3%] 
and a fractional rms deviation of 5%— 8%; the symmet- 
ric mass ratio 77 could be estimated with a bias within 
[-0.06, +0.04] and an rms deviation of 0.02-0.04; the 
chirp mass A4 could be estimated with a fractional bias 
within [—0.04%, +0.01%] and a fractional rms deviation 
of 0.7%-2%. We also proposed estimators for certain 



spin parameters of the double-spin system. For exam- 
ple, the parameter (S e g • Ln)/M 2 [where S c g is defined 
by Eq. J2EJ)], which is conserved when spin-spin effects 
can be neglected, could be estimated with a bias within 
[—0.01, +0.03] and an rms deviation of 0.07-0.13; the pa- 
rameter xtot = Stot/^f 2 could be estimated with a bias 
within [-0.005, +0.06] and an rms deviation of 0.10-0.20. 
However, since the mismatch metric has small compo- 
nents along the directions of these spin estimators, we 
expect that (at least for moderate S /N) statistical errors 
will always be dominant over the systematic errors dis- 
cussed here. We defer the study of statistical errors to a 
forthcoming paper |43j . 

In evaluating the performance of the quasi-physical 
single-spin template family, we have assumed a uniform 
distributions for the initial local parameters (spin and 
orientation angles) of the double-spin target model. It 
would be interesting in the future to redo our analyses 
assuming more realistic nonuniform distributions derived 
from astrophysical considerations. The only available re- 
sults for spin distributions in BH-BH binaries (unfortu- 
nately, with a single spin) and in NS-BH binaries were 
obtained using population-synthesis techniques 0> Hi| • 
For the case of binaries formed in globular clusters, there 
is no theoretical argument to suggest any particular spin 
distribution. 

Last, recent studies of PN spin-spin effects suggest 
that, for binaries with comparable masses, the two BH 
spins may have become roughly locked into a fixed rela- 
tive configuration by the time the GWs enter the band 
of good interferometer sensitivity. If these results are 
confirmed, they could provide preferred initial spin con- 
ditions, and by reducing the variability of GW signals, 
they could help to explain the good performance of our 
single-spin template family. We must however note that 
we were motivated in proposing our single-spin templates 
by the assumption that spin-spin effects never dominate, 
while for locking to occur spin-spin effects seem to be 
crucial |45|. 
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